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Abstract. The purpose of this paper is to describe geometrically discrete 
Lagrangian and Hamiltonian Mechanics on Lie groupoids. From a variational 
principle we derive the discrete Euler-Lagrange equations and we introduce a 
symplectic 2-section, which is preserved by the Lagrange evolution operator. 
In terms of the discrete Legendre transformations we define the Hamiltonian 
evolution operator which is a symplectic map with respect to the canonical 
symplectic 2-section on the prolongation of the dual of the Lie algebroid of 
the given groupoid. The equations we get include as particular cases the 
classical discrete Euler-Lagrange equations, the discrete Euler-Poincare and 
discrete Lagrange-Poincare equations. Our results can be important for the 
construction of geometric integrators for continuous Lagrangian systems. 
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1. Introduction 

During the last decade, much effort has been devoted to construction of geomet- 
ric integrators for Lagrangian systems using a discrete variational principle (see [21j 
and references therein). In particular, this effort has been concentrated for the case 
of discrete Lagrangian functions L on the cartesian product Q x Q of a differentiable 
manifold. This cartesian product plays the role of a "discretized version" of the 
standard velocity phase space TQ. Applying a natural discrete variational princi- 
ple, one obtains a second order recursion operator ^ : Q x Q — > Q x Q assigning 
to each input pair {x,y) the output pair {y,z). When the discrete Lagrangian is 
an approximation of a continuous Lagrangian function (more appropriately, when 
the discrete Lagrangian approximates the integral action for L) we obtain a nu- 
merical integrator which inherits some of the geometric properties of the contin- 
uous Lagrangian (symplecticity, momentum preservation). Although this type of 
geometric integrators have been mainly considered for conservative systems, the 
extension to geometric integrators for more involved situations is relatively easy, 
since, in some sense, many of the constructions mimic the corresponding ones for 
the continuous counterpart. In this sense, it has been recently shown how dis- 
crete variational mechanics can include forced or dissipative systems, holonomic 
constraints, explicitely time-dependent systems, frictional contact, nonholonomic 
constraints, multisymplectic fields theories... All these geometric integrators have 
demonstrated, in worked examples, an exceptionally good longtime behavior and 
obviously this research is of great interest for numerical and geometric considera- 
tions (see |H|). 

On the other hand, Moser and Veselov consider also discrete Lagrangian 
systems evolving on a Lie group. All this examples leads to A. Weinstein [31] to 
study discrete mechanics on Lie groupoids, which is a structure that includes as 
particular examples the case of cartesian products Q x Q as well as Lie groups. 

A Lie groupoid G is a natural generalization of the concept of a Lie group, 
where now not all elements are composable. The product gig2 of two elements is 
only defined on the set of composable pairs G2 — {(g, h) d G x G \ (3{g) = a(/i)} 
where a : G — > M and /3 : G — > M are the source and target maps over a 
base manifold M. This concept was introduced in differential geometry by Ch. 
Ereshmann in the 1950's. The infinitesimal version of a Lie groupoid G is the Lie 
algebroid AG — > M, which is the restriction of the vertical bundle of a to the 
submanifold of the identities. 

We may thought a Lie algebroid A over a manifold M, with projection r : 
A M, as a generalized version of the tangent bundle to M. The geometry and 
dynamics on Lie algebroids have been extensively studied during the past years. 
In particular, one of the authors of this paper (see [12]) developed a geometric 
formalism of mechanics on Lie algebroids similar to Klein's formalism [11] of the 
ordinary Lagrangian mechanics and more recently a description of the Hamiltonian 
dynamics on a Lie algebroid was given in [THj^ (see also [IB])- 

The key concept in this theory is the prolongation, V^A, of the Lie algebroid over 
the fibred projection r (for the Lagrangian formalism) and the prolongation, A, 
over the dual fibred projection t* : A* — > M (for the Hamiltonian formalism). See 
[l4] for more details. Of course, when the Lie algebroid is A = TQ we obtain that 
VA = T{TQ) and A = T{T*Q), recovering the classical case. An alternative 
approach, using the linear Poisson structure on A* and the canonical isomorphism 
between T*A and T*A* was discussed in [7]. 

Taking as starting point the results by A. Weinstein [31], we elucidate in this 
paper the geometry of Lagrangian systems on Lie groupoids and its Hamiltonian 
counterpart. Weinstein gave a variational derivation of the discrete Euler-Lagrange 
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equations for a Lagrangian L : G ^ K on a Lie groupoid G. We show that the 
appropriate space to develop a geometric formahsm for these equations is the Lie 
algebroid V^G = VP 0g Vet — > G (see section 3 for the definition of the Lie 
algebroid structure). Note that V^G is the total space of the prolongation of the 
Lie groupoid G over the vector bundle projection r : AG M, and that the Lie 
algebroid of V^G is just the prolongation V'^{AG) (the space were the continuous 
Lagrangian Mechanics is developed). Using the Lie algebroid structure of V^G we 
may describe discrete Mechanics on the Lie groupoid G. In particular, 

- We give a variational derivation of the discrete Euler-Lagrange equations: 

X{g){L)-X{h){L) = Q 

for every section X of AG, where the right or left arrow denotes the induced 
right and left-invariant vector field on G. 

- We introduce two Poincare-Cartan 1-sections and O^, and an unique 
Poincare-Cartan 2-section, fii, on the Lie algebroid P'^G — > G. 

- We study the discrete Lagrangian evolution operator ^ : G — > G and its 
preservation properties. In particular, we prove that {V^ S^^S^Y^l = ^l, 
where V^S, is the natural prolongation of ^ to V^G. 

- Reduction theory is stablished in terms of morphisms of Lie groupoids. 

- The associated Hamiltonian formalism is developed using the discrete Le- 
gendre transformations ¥+ L : G A*G and ¥^ L : G A*G. 

- A complete characterization of the regularity of a Lagrangian on a Lie 
groupoid is given in terms of the symplecticity of Ql or, alternatively, 
in terms of the regularity of the discrete Legendre transformations. In 
particular, Theorem l4. 1 31 solves the question posed by Weinstein [ST about 
the regularity conditions for a discrete Lagrangian function on more general 
Lie groupoids than the cartesian product Q x Q. In the regular case, we 
define the Hamiltonian evolution operator and we prove that it defines a 
symplectic map. 

- We prove a Noether's theorem for discrete Mechanics on Lie groupoids. 

- Finally, some illustrative examples are shown, for instance, discrete Me- 
chanics on the cartesian product QxQ, on Lie groups (discrete Lie-Poisson 
equations), on action Lie groupoids (discrete Euler-Poincare equations) 
and on gauge or Atiyah Lie groupoids (discrete Lagrange-Poincare equa- 
tions). 

We expect that the results of this paper could be relevant in the construction of 
new geometric integrators, in particular, for the numerical integration of dynamical 
systems with symmetry. 

The paper is structured as follows. In Section 2 we review some basic results on 
Lie algebroids and Lie groupoids. Section 3 is devoted to study the Lie algebroid 
structure of the vector bundle V^G = VP (Bg Va G. The main results of the 
paper appear in Section 4, where the geometric structure of discrete Mechanics 
on Lie groupoids is given. Finally, in Section 5, we study several examples of the 
theory. 

2. Lie algebroids and Lie groupoids 

2.1. Lie algebroids. A Lie algebroid A over a manifold M is a real vector bundle 
T : A ^ M together with a Lie bracket |-, •] on the space T{t) of the global cross 
sections of r : ^ ^ M and a bundle map p : A ^ TM, called the anchor map, 
such that if we also denote by p : T{t) X(M) the homomorphism of C°°(A/)- 
modules induced by the anchor map then 

lXJYj = flX,Yj+piX)if)Y, (2.1) 
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for X, r e r(T) and / e C°°(M) (see [TT]). 

If X,Y,Z e T{t) and / G C°°{M) then, using and the fact that |-, •] is a 
Lie bracket, we obtain that 

IIX, Y], fZ] = filX, lY, Z]] - lY, IX, Z]]) + [p{X),p{Y)]{f)Z. (2.2) 

On the other hand, from (|2.ip . it follows that 

IIX, Yl fZ] fllX, YlZj+ pfX, Y]if)Z. (2.3) 

Thus, using (|2.2p . (|2.3p and the fact that |-, •] is a Lie bracket, we conclude that 

plX,Yj - [p{X),p{Y)], 

that is, p : r(r) X{M) is a homomorphism between the Lie algebras (r(T), |-, •]) 
and (X(M), [•,•]). 

If {A, |-, •], p) is a Lie algebroid over M, one may define the differential of A, 

d : r(A'=r*) ^ r(A'=+V*), as follows 



k 



dpixo, . . . , Xfe) = }_^{-iypix,){fi{Xo, . . . , x„ . . . , Xfe)) 

1=0 (2.4) 
+ Xj,Xo, . . . ,X;, . . . , X;-, . . . ,Xfe), 

■i<i 

for /i e r(A'^T*) and Xq, . . . , Xfe G r(r). c? is a cohomology operator, that is, — 0. 
In particular, if / : Af — > M is a real smooth function then df{X) — p{X)f, for 
X G r(r). We may also define the Lie derivative with respect to a section X of 
A as the operator Cx : r(A''A*) — > r(A''^*) given by £x ^ ix o d + d o ix (for 
more details, see [Vf]). 

Trivial examples of Lie algebroids are a real Lie algebra q of finite dimension (in 
this case, the base space is a single point) and the tangent bundle TM of a manifold 
M. Other examples of Lie algebroids are: i) the vertical bundle (tp)|v^ : Vn P 
of a fibration it : P ^ M (and, in general, the tangent vectors to a foliation of finite 
dimension on a manifold P); ii) the Atiyah algebroid associated with a prin- 
cipal G-bundle (see [HJIT^); iii) the prolongation V^A of a Lie algebroid 
A over a fibration tt : P ^ M (see [9l[T4j) and iv) the action Lie algebroid 
AtK f over a map f : M' ^ M (see [Ti]V 

Now, let {A, |-, ■], p) (resp., {A', |-, •]', p')) be a Lie algebroid over a manifold M 
(resp., M') and suppose that ^> : A ^ A' is a vector bundle morphism over the map 
^'o : M — > M'. Then, the pair (^f, ^fo) is said to be a Lie algebroid morphism if 

(i((«', *o)*0') = (*, *o)*(rfV'), for aU 0' G r{A''{A')*) and for aU k, (2.5) 

where d (resp., d') is the differential of the Lie algebroid A (resp.. A') (see !14 ). In 
the particular case when M — M' and = Id then (|2.5p holds if and only if 

Ivi'oX,vi/or]' = vi/|x,ri, p'{^x)^p{X), iorX,YeT{T). 

2.2. Lie groupoids. In this Section, we will recall the definition of a Lie groupoid 
and some generalities about them are explained (for more details, see [51 117|). 
A groupoid over a set M is a set G together with the following structural maps: 
• A pair of maps a : G — » Af, the source, and j3 : G ^ M, the target. 

Thus, an element 5 G G is thought as an arrow from x = a{g) to y = (3{g) 
in M 



X = 0,(9) = 0(g) 

The maps a and /? define the set of composable pairs 



G2 = {{g,h)eGxG/l3{g)^a{h)}. 
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• A multiplication m : G2 ^ G, to be denoted simply by m{g, h) — gh, 
such that 

— a{gh) — a{g) and f3{gh) — (3{h). 

- g{hk) = {gh)k. 

If g is an arrow from x = a{g) to y = f3{g) — a{h) and h is an arrow from 
y to z = (3{h) then gh is the composite arrow from a; to 2; 




x^a(g)^a(gh) y ^ 13(g) = a(h) z = I3{h) = (3{gh) 

• An identity section e : M G such that 

- e(a(5)).9 = 9 and ge{(3{g)) = g. 

• An inversion map i : G ^ G, tohe denoted simply by i{g) = g~^, such 
that 

- 9^^9 = and gg-^ =_e{a{g)). 

9 



0(9) = a(g-^) 




A groupoid G over a set M will be denoted simply by the symbol G =t M. 

The groupoid G =t M is said to be a Lie groupoid if G and AI are manifolds and 
all the structural maps are differentiable with a and /3 differentiable submersions. 
If G ^ M is a Lie groupoid then m is a submersion, e is an immersion and i is 
a diffeomorphism. Moreover, if x £ M, a^-^{x) (resp., /3~^(a;)) will be said the 
a-fiber (resp., the (3-fiber) of x. 

On the other hand, if g G G then the left-translation by g £ G and the 
right-translation by g are the diffeomorphisms 

;g:a"i(/3(g))^a-i(a(g)) ; h 



l,{h) : 



gh, 

hg. 



Note that 1^^ 



Ig-i and Tg ^ 



A vector field X on G is said to be left-invariant (resp., right-invariant) if 
it is tangent to the fibers of a (resp., /?) and X[gh) = {Thlg){Xh) (resp., X[gh) = 
(V^)(X(3))), for {g,h)&G2. 

Now, we will recall the definition of the Lie algebroid associated with G. 

We consider the vector bundle r : AG M, whose fiber at a point x G M 
is Ax:G = V^(^x)0^ = -?^er(Tg(3,')a). It is easy to prove that there exists a bijection 
between the space T{t) and the set of left-invariant (resp., right-invariant) vector 
fields on G. If X is a section of r : AG — > M, the corresponding left-invariant 
(resp., right-invariant) vector field on G will be denoted X (resp., X), where 



X{g) = (T,(^(g))Z,)(X(/3(.g))), 

^(.9) = ~{Teia(g))rg){{T,^a(g))i){X{a{g)))), 



(2.6) 
(2.7) 



for g G G. Using the above facts, we may introduce a Lie algebroid structure 
(|-, -JjP) on AG, which is defined by 



[X, Yj = [X, Y], p{X){x) = (T,(,)/3)(X(x)), 
for X,Y e V{t) and x e M. Note that 



IX, Yj 

TioX^ 
(for more details, see [il fT7]). 



---[X,Y], [X,Y]^0, 
-X o i, Ti o X — —X o i, 



(2.8) 

(2.9) 
(2.10) 
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Given two Lie groupoids G ^ M and G" ^ M' , a morphism of Lie groupoids 
is a smooth map ^ : G ^ G' such that 

{g,h) e G2 =^ m9),<i>{h)) e {G'h 

and 

A morphism of Lie groupoids $ : G ^ G" induces a smooth map : M ^ M' in 
such a way that 

a, (3 and e (resp., a', /?' and e') being the source, the target and the identity section 
of G (resp., G'). 

Suppose that ($, $0) is a morphism between the Lie groupoids G =t M and 
G' =4 A/' and that t : AG ^ M (resp., r' : AG' M') is the Lie algebroid of 
G (resp., G'). Then, if x G M we may consider the hnear map Ax{^) : A^G 
^■s>oio!:)G' defined by 

AA^){ve{cc)) = (r,(^)$)(i;,(^)), for e ^a.G. (2.11) 

In fact, we have that the pair (A($), $0) is a morphism between the Lie algebroids 
T : AG ^ M and t' : AG' ^ M' (see (T^). 

Next, we will present some examples of Lie groupoids. 

1. - Lie groups. Any Lie group G is a Lie groupoid over {e}, the identity element 
of G. The Lie algebroid associated with G is just the Lie algebra g of G. 

2. - The pair or banal groupoid. Let AI be a manifold. The product manifold 
M X M is a Lie groupoid over M in the following way: a is the projection onto 
the first factor and /3 is the projection onto the second factor; e(x) — {x, x), for all 
X G M, m{{x, y), {y, z)) = {x, z), for {x, y),{y,z) € M x M and i{x, y) = (y, x). 
M X M ^ M is called the pair or banal groupoid. If a: is a point of M, it follows 
that 

= {Ox} X T,M C TxM x T,M = T(,,,)(Af x M). 
Thus, the linear maps 

*x : T^M Ve(x)a, (Ox, fx), 

induce an isomorphism (over the identity of M) between the Lie algebroids Tm : 
TM ^ M and T : A{M x M) ^ M. 

3. - The Lie groupoid associated with a fibration. Let tt : P ^ A/ be a 

fibration, that is, tt is a surjectivc submersion and denote by Gtt the subset of P xP 
given by 

G^ ^ {{p,p') e P X Pl-Kip) = ^{p')]. 

Then, Gtt is a Lie groupoid over P and the structural maps a^r, P-m "ni-m e-tt and 
are the restrictions to Gtt of the structural maps of the pair groupoid P x P ^ P. 
If p is a point of P it follows that 

V^Ap)"- = {(Op,i"p) e TpP X TpP/{Tpn){Yp) = 0}. 

Thus, if (Tp)|y^ : Vt: P is the vertical bundle to vr then the linear maps 

{■f^)p : VpTT V, „(p)Q!Trj Yp > {Op,Yp) 

induce an isomorphism (over the identity of A/) between the Lie algebroids {tp)^y^ : 
Ftt ^ P and r : AG^ P. 

4. - Atiyah or gauge groupoids. Let p : Q — * A/ be a principal G-bundle. 
Then, the free action, <& : G x Q — > Q, {g,q) ^ $(5, q) = gq, of G on Q induces, in 
a natural way, a free action <I>x$:Gx(Qx(5)^(5x(3ofGon(3x(5 given 
by ($ X $)(g, (g, q')) ~ {gq, gq'), for g G G and {q, q') e Q x Q. Moreover, one may 
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consider the quotient manifold {Q x Q)/G and it admits a Lie groupoid structure 
over M with structural maps given by 



a : (Q X Q)/G — > M 

/3 : (Q X Q)/G ^ M 

e:M^{QxQ)/G 

m:{{Qx Q)/G)2 ^ (Q x Q)/G 

~i:{Qx Q)/G -^{Qx Q)/G 



([(?, <?')], [(V, 9")])- 

[{q,q')]~^[{q\q)]. 



X 

[(m,q'% 



This Lie groupoid is called the Atiyah (gauge) groupoid associated with the 
principal G-bundle p : Q ^ M (see [16]). 

If a; is a point of M such that p{q) = x, with q E Q, and pqxq '■ Q x Q ^ 
{Q X Q)/G is the canonical projection then it is clear that 



(T^(?.g)PQxQ)({Ojxr,Q). 



. -p^Q ^ p 

V'G 



m 



: P 

: [V''G)2 — > V'G 
VG — > VG 



Thus, if tq I G : TQ/G ^ M is the Atiyah algebroid associated with the principal 
G-bundlc p : G ^ M then the linear maps 

(TQ/G)^ Vi(^)a ; [vq\ {Tiq^q)PQxQ){Oq,Vq), with Vg e TgQ, 

induce an isomorphism (over the identity of M) between the Lie algebroids r : 
MQ X Q)/G) M and tq\G : TQ/G M. 

5. - The prolongation of a Lie groupoid over a fibration. Given a Lie 
groupoid G ^ M and a fibration tt : P ^ M, we consider the set 

P^G = P^x„G px^ P = {{p,g,p') ePxGx P/nip) = a{g), I3{g) = tt{p')}. 

Then, V^G is a Lie groupoid over P with structural maps given by 

a-:V^G^P ■ ip,g,p')^p, 

ip,9,p') — 'P', 
P — ' (P,e(7r(p)),p), 
{(j>,g,p'),{p' ,h,p")) ^ {p,gh,p"), 

{p,g,p') — ' {p',9^^,p)- 

V^G is called the prolongation of G over tt : P M. 

Now, denote by r : AG M the Lie algebroid of G, by A{V^G) the Lie algebroid 
of V^G and by V'^{AG) the prolongation of r : AG M over the fibration tt. If 
p € P and m = tt{p), then it follows that 

ApiP^G) = {(Op, Xp) e TpP X A„,G x Tj,P/ {Tj,n){Xp) = (r,(„)/3)(z;,(„))} 

and, thus, one may consider the linear isomorphism 

{■^nj,:Aj,{V^G)^r;[AG), (Op, i;,(„,), Xp) ^ («,(„,), Xp). (2.12) 

In addition, one may prove that the maps (^''^)p, p E P, induce an isomorphism 
-^"^ : AiV'G) ViAG) between the Lie algebroids AiV'G) and V^iAG) (for 
more details, see 

6. - Action Lie groupoids. Let G =t M be a Lie groupoid and tt : P ^ M he 
a smooth map. If P ^Xq, G — {(p, g) G P x G/tt{p) — a{g)} then a right action of 
G on TT is a smooth map 

P^XaG^P, {p,g)^pg, 

which satisfies the following relations 

■^{pg) = Pig), ioT {p,g) € P ^Xc,G, 
{pg)h = p{gh), for {p,g) e P ^x^ G and [g, h) e G2, and 
pt{TT{p)) = p, for p e P. 
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Given such an action one constructs the action Lie groupoid P ^ x „ G over 
P by defining 



: P G — > P 

Pit '■ P -K^ a G > P 

Ctt • P * P TT^ a G 

rhr, : {P ^Xa G)2 — > P ^XaG 

i-n '• P TT^a G > P -j^X ^ G 



(P-, 9) — ^ P9-, 

{{P,9), iP9,h)) — > {p,gh), 
iP,9) — ' iP9,9^^)- 
Now, if p e P, we consider the map p ■ : a~^(7r(p)) P given by 

P- {9)= P9- 

Then, ii t : AG M is the Lie algebroid of G, the R-linear map $ : T{t) X{P) 
defined by 

$(X)(p) = (r,(,(p))p ■){X{7r{p))), for X e r(T) and p G P, 

induces an action of AG on vr : P ^ M. In addition, the Lie algebroid associated 
with the Lie groupoid PttXqG^P is the action Lie algebroid AG k tt (for more 
details, see j9]). 

3. Lie algebroid structure on the vector bundle vr"^ : V^G G 

Let G ^ M be a Lie groupoid with structural maps 

a,(3:G^ M, e : M ^ G, i:G^G, m : G2 ^ G. 

Suppose that r : AG — ^ M is the Lie algebroid of G and that V^G is the prolon- 
gation of G over the fibration r : AG M (see Example 5 in Section 12. 2p , that 
is, 

P^G = ylG^x<,G^x^ylG. 
V^G is a Lie groupoid over AG and we may define the bijective map : V^G 
V(3 ©G as follows 

for (Me(Q(g)),5','ye(/3(g))) S ^a(g)G X G X ^/3(g)G. Thus, VP^aVa is a Lie groupoid 
over AG (this Lie groupoid was considered by Saunders [3^). We remark that the 
Lie algebroid of V^G = V(3 ®g Vet =^ AG is isomorphic to the prolongation of 
AG over r : AG M and that the prolongation of a Lie algebroid A over the 
vector bundle projection t : A —>■ M plays an important role in the description of 
Lagrangian Mechanics on A (see [HI 122) ). 

On the other hand, note that V^G = Vp ©g is a real vector bundle over 
G. In this section, we will prove that the vector bundle tt'^ : V^G G admits an 
integrable Lie algebroid structure. In other words, we will prove that there exists 
a Lie groupoid =4 G over G such that the Lie algebroid AH is isomorphic to the 
real vector bundle tt"^ : V^G — > G. In addition, we will see that the Lie groupoid 
H is isomorphic to the prolongations of G over a and p. 

It is clear that the Lie algebroids of the Lie groupoids over G 

Gp = {(.g, h)eGx G/p{g) = p{h)}, G„ = {(r, s) € G x G/a(r) = a(s)}, 

are just VP G and Va — > G, respectively. This fact suggests to consider the 
following manifold 

Gp-kGa = {{{g, h), (r, s)) £Gp X Ga/Ppig, h) = aa{r, s)}, 

where Pp : Gp ^ G (respectively, aa ■ Ga ^ G) is the target (respectively, the 
source) of the Lie groupoid Gp ^ G (respectively, Gq ^ G) . 
We will identify the space Gp ★ Gq with 

{(5, /i, s) e G X G X G/P{g) = P{h), a{h) = a{s)}. 
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ig,h,s) — > g, 
{g,h,s) — > s, 

9 * {9: 9: 9): (3-1) 

{{g,h,s),{s,h',s'))^^{g,h's-^h, s'), 
{g,h,s) — {s,gh-'^s,g). 



This space admits a Lie groupoid structure over G with structural maps given by 
a/3a Gp-kGa — > G 

Ppa ■ Gp -k Ga > G 

epa ■ G > G(3 ★ Ga 

m/Ba ■ {Gp -k Ga)2 * Gp ★ Gq 

ipa ■ Gp-kGa > Gp -k Ga 

Note that 

jp:Gp — >GpkGa ; {g.h) — > jp{g,h) = {g,h,h), 

ja-Ga > Gp-kGa ] {h, s) > j a{h, s) = {h, h, s) , 

are Lie groupoid morphisms and that the map 

mpa{jp,ja) ■■ Gp-kGa^ GpkcGa', {g,h,s) mpa{j p {g , h) , ja{h, s)) 

is just the identity map. This imphes that {Gp,Ga) is a matched pair of Lie 
groupoids in the sense of Mackenzie [18] (see also [25] ) . 

Denote by (|-, -Jjp) the Lie algebroid structure on r : AG M. 

Theorem 3.1. Let A{Gp -k Ga) ^ G be the Lie algebroid of the Lie groupoid 
Gpi^Ga^G. Then: 

(i) The vector bundles A{Gp k G a) ^ G and n'' : VG '^VI3®GVa^ G are 
isomorphic. Thus, the vector bundle tt'^ : V^G = VP (Bg ^ol — > G admits 
a Lie algebroid structure. 

(ii) The anchor map o/ tt"^ : V^G ^ V (3 ®g —> G is given by 

P^^''{Xg,Yg) =Xg+ Yg , for {X g , Yg) £ Vg ^ ffi Vg U , (3.2) 

and the Lie bracket |-,-]^ '"^ on the space r(7r'^) is characterized by the 
following relation 

I(X, Y), {X', r)f^^ = i-lX, X'llYX]), (3.3) 
for X,Y,X',Y' eVir). 

Proof, (i) If 5 e G then, from p.ip . we deduce that the vector space V^^^i^g-^apa 
may be described as follows 

Ve,^ig)apa = {(O,,,^,,, Zg) G TgG X TgG X TgG/Xg e VgP, (T<,«)(Xg) = iTga)iZg)} 

- {{Xg, Zg) e TgG X TgG / X g £ Vg (3 , (r,«)(X,) = iTga)iZg)}. 

Now, we will define the linear map 'i'g : Vti^,^(g)Oipa Vgf3 Vga = 'PgG by 

^g{Xg,Zg) = {Xg,Zg-Xg). (3.4) 

It is clear that ^'g is a linear isomorphism and 

%\Xg,Yg) = iXg,Xg + Yg) , fOT [X g , Yg) £ Vg ^ © VgU ^ G . (3.5) 

Therefore, the collection of the maps g, g E G , induces a vector bundle isomor- 
phism : A{Gp ★ Ga) VG = VP ®g Va over the identity of G. 

(ii) A direct computation, using ()3.ip . proves that the linear map T^fj^(g)f3pa '■ 
^e3=(9)"/3" ^ TgG is given by 

iTe^^^ig)Ppa){Xg,Zg) = Zg. (3.6) 

Consequently, from (|2.8p . p.5p and p.6p . we deduce that (|3.2p holds. 
Next, we will prove (j3.3p . 
Using (|3.5p . it follows that 

(vp-i o (X, r ))(g) = (0„ X{g), X{g) + Y{g)) = (X{g), X{g) + Y{g)), 
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for g E G. Denote by^* ^ o [X ,Y) the corresponding left-invariant vector field on 
G/3 * Ga- Then, from (|2.6p and (|3.ip . we have that 



o {X, Y)ig, h, s) = (Og, X{h),X{s) + r (s)), for {g, h, s) £ G/3 * G„. 
Thus, using (j2.8p and ()2.9p . we conclude that 



Therefore, we obtain that ()3.3|) holds. □ 
The above diagram shows the Lie groupoid and Lie algebroid structures of V^G: 



r^G 




=^AG 




Tf3 








T 





TM 



^ M 



Given a section X of AG — > A/, we define the sections X^^^^\ (the /3 

and a- lifts) and X^^'^'^ (the complete hft) of X to tt^ : VG — > G as follows: 

X(i'°)(5) = (X(g),0,), X(«'i)(g) = (0„X(g)) and X^'^^\g) = {-X{g),X[g)) 

We can easily see that 

|jf(0,l) y(0,l)jP"G ^ IX ^ ~ ^ 

and, as a consequence, 

1^(1,1) y(0,l)j-p-G ^ yj(0,l) a^^d ^ (3.8j 

Remark 3.2. From Theorem 13. 11 we deduce that the canonical inclusions 

{Id, 0) : y/3 ^ V^G ^ VP ©g Va, (0, /rf) : Va VG ^ Vf3 ®g Va, 

are Lie algebroid morphisms over the identity of G. In other words, {VP, Va) is 
a matched pair of Lie algebroids in the sense of Mokri [21] . This fact directly 
follows using the following general theorem (see [25 ): if {G,H) is a matched pair 
of Lie groupoids then {AG, AH) is a matched pair of Lie algebroids. o 

Next, we will consider the prolongation V'^G of the Lie groupoid G over the 
target p-.G^M. We recaU that 

V^G ^Gpx^G fixp G = {{g, h,s)eGxGx G/P{g) = a{h), P{h) = P{s)}, 

and that V^G is a Lie groupoid over G with structural maps 



: pf^G — > G 
pP : V^G — > G 
ef^ -.G — > V/^G 



{V^G), 



V^G 



iP : rl^G — y V'^G 



(3.9) 



{g,h,s) — > g, 
{g,h,s) — > s, 

9 — ' i9,<^iP{g)),g), 

{{g, h, s), {s, t, u)) — > {g, ht, u), 
{g,h,s) — > {s,h~^,g). 

Moreover, we also have that the Lie algebroid of V^G may be identified with the 
prolongation V^{AG) of AG over P : G — > M . We remark that 

V^{AG) = {(i;,(^(,)),X,) e A^(,)G X T,G/(T,(^(,))/3)(z;,(^(,))) = (r,/?)(X,)} 

for g eG. 
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Theorem 3.3. Let <^f^ : G p -k P^G be the map defined by 

= (3.10) 

for [g, h, s) G Gp -k Ga- Then: 

(i) is a Lie groupoid isomorphism over the identity of G. 
ill) If A{<^f^) : A{Gp ★ Go) — > A[V^G) is the corresponding Lie algebroid 
isomorphism then, under the identifications 

A{Gp * Ga) = V^G ^ V(3 ©G Va, AiV/^G) ^ V^{AG), 

A{^^) is given by 

Agi^f'){Xg,Yg) = {{T ,1 .) (Y ,) , X , + Y ,) , (3.11) 

for {Xg^Yg) G Vgl3 ® Vgtt , whcre Ig-i : a^^{a{g)) a~^{f3{g)) is the 
left-translation by g~^ . 

Proof, (i) A direct computation, using (|3.ip and (|3.9p . proves the result, 
(ii) If 5 G G we have that 

AgiGp * G„) = V,^^^g)af3a = {{Og, Xg, Zg) e TgG X TgG X TgG/Xg G Vg!3, 

{Tga){Xg) = {Tga){Zg)}, 

AgipPG) = V,,(g)aP = {{%,v,(pig)),Yg) G TgG X Ap(g)G x TgG/ 

{TgP){Yg) = {T,(P(g))l3){v,^P(g)))}. 

Now, if {Qg,Xg,Zg) G V,^^ .{g)Ci^l3a then, from (|3.10p . we deduce that 
iTe,^ig)'i'^)iOg,Xg,Zg) = (r,^,„(,)$^)(0„Og,Z, -X,) + (T,^,^(,)Ci>^)(Og,X„X,) 

= {Og, {Tglg-.){Zg - X g) , Zg - X g) + ( „ ( , ) <i> ) ( Q,, X„ X, ) . 

On the other hand, suppose that (3{g) — x ^ M and that 7 : (—£,£) /3~^(x) is a 
curve in (3~'-^{x) such that 7(0) = g and 7'(0) = Xg. Then, one may consider the 
curve 7 : (— e, e) Gfj ★ Ga on Gp ★ Ga given by 

l(.t) = {g,l{t),j{t)) 

and it foUows that 

7(0) = (5,5,5), j'{0) = iOg,Xg,Xg). 

Moreover, we obtain that 

7(t) = i^Poj)it) = ig,e{x),j{t)), for all t 

and thus 

7'(0) = (0g, 0,(0(3)), X,). 

This proves that 

{T,^^^g^<^>^){Og,Xg,Zg) = {O g , {Tgl g - ^ ) {Z g - X g) , Z g) . (3.12) 

Finally, using (PH]) . ((TT^ . (^3)) and ((XT^ . we deduce that (PTT|l holds. □ 

Next, we will consider the prolongation P^G of the Lie groupoid G over the 
source a : G — > M. We recall that 

V^G = G aX-cG pXaG^ {(5, /i, s) G G x G X G/a(5) a{h),P{h) = a(s)} 

and that P^G is a Lie groupoid over G with structural maps 

a^-.V^G—^G ; {g,h,s)^g, 

f}°':P°'G — >G ■ {g,h,s) — > s, 

,^:G~^P''G ; g ^ (g , eiaig)) , g) , 

m'- -.{P'-Gh^P'^G ; {{g,h,s),{s,t,u)) ^ {g,ht,u), 

-.P^G ^P'^G ■ {g,h,s)~^{s,h~\g). 
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Moreover, we also have that the Lie algebroid of V^G may be identified with 
the prolongation V°'{AG) of AG over a : G ^ M. We remark that 

for g e G. 

Theorem 3.4. Lei : Gfj -k G^ ^ V^G be the map defined by 

$"(5,/i,s) = (5r,5r/i-\s), 

/or {g, h, s) G G/3 ★ Gq,. Then: 

(i) is a Lie groupoid isomorphism over the identity of G. 

(ii) // : A(G/3 ★ Gq) ^ AiV^G) is the corresponding Lie algebroid 
isomorphism then, under the canonical identifications 

A{Gf} * Ga) = V^G ^ Vf3 ®G Va, A(7'"G) =^ ViAG), 

is given by 

Ag{^^){Xg,Yg) = {Tg {l O r - ^ ) {X g) , X g + Yg) , (3.13) 

for {Xg,Yg) G K,/3©K,a, w/iere r^-i : /3-i(/3(g)) /3"Ha(5)) «s 
right-translation by g^^ . 

Proof. Proceeding as in the proof of Theorem 13.31 we deduce the result. □ 

4. Mechanics on Lie Groupoids 

In this section, we introduce Lagrangian (Hamiltonian) Mechanics on an ar- 
bitrary Lie groupoid and we will also analyze its geometrical properties. This 
construction may be considered as a discrete version of the construction of the 
Lagrangian (Hamiltonian) Mechanics on Lie algebroids proposed in [52] (see also 
[131 [53]). We first discuss discrete Euler-Lagrange equations following a similar ap- 
proach to [3T], using a variational procedure. Secondly, we intrinsically define and 
discuss the discrete Poincare-Cartan sections, Legendre transformations, regularity 
of the Lagrangian and Noether's theorem. 

4.1. Discrete Euler-Lagrange equations. Let G be a Lie groupoid with struc- 
tural maps 

a,(3:G-^M, e: M -^G, i:G^G, m : G2 ^ G. 

Denote by r : AG M the Lie algebroid of G. 

A discrete Lagrangian is a function L : G — > R. Fixed g E G, we define the 
set of admissible sequences with values in G: 

Cg = {(.91, . . . , 5Jv) e G^ / (gfc, .gfc+i) e G2 for fc = 1, . . . , iV- 1 and gi . . . .g„ = g}. 

Given a tangent vector at {gi, . . . ,gN) to the manifold , we may write it as 
the tangent vector at t = of a curve in , t G (— e, e) C R — > c{t) which passes 
through (51, ... , gpf) a.t t — 0. This type of curves is of the form 

c{t) = (gihi{t),h^^{t)g2h2{t), h^\^{t)gN ^ihN -i{t), h'^\^{t)gN) 

where hk{t) £ a~'^{(3{gk)), for aU t, and /ifc(O) = e(/3(5fc)) for fc = 1, . . . , iV 1. 
Therefore, we may identify the tangent space to at (51, . . . , g^) with 

r(g,,...,g„)Cf = {(wi, . . .,VN^i) I vk e A^^G and Xk = /3(.9fe), 1 < fc < - 1} . 

Observe that each Vk is the tangent vector to the a-vertical curve hk at t = 0. 

The curve c is called a variation of {gi, . . . , gj^) and (vi, V2, ■ ■ ■ , vn-i) is called 
an infinitesimal variation of {gi, . . . , g^). 
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Define the discrete action sum associated to the discrete Lagrangian L : G — > 

M 

SL : — > R 

N 
k=l 

We now proceed, as in the continuous case, to derive the discrete equations 
of motion applying Hamilton's principle of critical action. For it, we consider 
variations of the discrete action sum. 



Definition 4.1 (Discrete Hamilton's principle 31J). Given g £ G, an admissible 



sequence {gi, . . . ,gN) S is a solution of the Lagrangian system determined by 



L : G > R if and only if {gi, . . . , gjf) is a critical point of SL. 

Fist of all, in order to characterize the critical points, we need to calculate: 



d 
di 



d 



t=o dt 



5L(c(i)) - - {Ligihiit)) + Lih^\t)g2h2{t)) 



+ ... + L{hj^\^{t)gN-ihN-i{t)) + L{h-\{t)gN)) • 

Therefore, 

N-l 

SL(r(f^] = 

dt 



SL{c{t)) = (d°(i o lg^){e{xk)){vk) + d°(L o r,,^, o z){e{xu)){vu)) 

* k=l 



where d° is the standard differential on G, i.e., d° is the differential of the Lie 

algebroid tq '■ TG G. Since the critical condition is — 

___ dt t=o 

applying (|2.6p and (|2.7p . we may rewrite this condition as 



SL{c{t)) = then. 



Af-l _ _ N-1 

= 

k=l k=l 



jv — 1 jv — 1 

J2 [Xki9km-X,{g,+,){L)] = J2 [{dL,xi'-'^){g,) ^ {dL,xi'''%,+,) 



where d is the differential of the Lie algebroid tt'^ : V^G = V(3 ©g Vet — > G and 
Xk is a section of r : AG M such that Xk{xk) ~ Vk- 

For iV = 2 we obtain that (51, 32) € G2 (with (3{gi) — a(.92) = x) is a solution if 

d° [Lolg^+Lo rg^ o i] {e{x))\A^G = 

or, alternatively, 

X{g,){L)-X{g2){L)=(} 

for every section X of AG. These equations will be called discrete Euler- Lagrange 
equations. 

Thus, we may define the discrete Euler-Lagrange operator. 

i'DEL^ : G2 ^ A*G , 

where A*G is the dual of AG. This operator is given by 

^DEL-^(5, h) ^ d^^ [L o Ig + L o rh o i] {e{x))\A^G 

with I3{g) = a{h) — x. 

In conclusion, we have characterized the solutions of the Lagrangian system 
determined by L : G — > M as the sequences (gi, . . . ,gN), with {gk,gk+i) G G2, for 
each fce{l,...,Af — 1}, and 

-DDELi(.9fc,5fe+i) = 0, 1 < A: < iV - 1. 
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4.2. Discrete Poincare-Cartan sections. Given a Lagrangian function L : G — : 
K, we will study the geometrical properties of the discrete Euler-Lagrange equa- 
tions. 

Consider the Lie algebroid tt'^ : P^G = V(3 (Bg Vol — > G, and define the 
Poincare-Cartan 1-sections 0^,8j G r((7r'^)*) as follows 

Ql(9){X,,Y,)^-X,{L), e+{g){Xg,Yg)^Yg{L), (4.1) 

for each g ^ G and {Xg, Yg) G VgP © Vga. From the definition, we have that 

eiig){X('-'H9))-^X{g)iL) and e^ig)iX(°''H9)) ^ 0, 

and similarly 

e+(5)(X("-i)(5)) = X(5)(i) and e+ig){X^'^°\g))^0, 
for X e T{t). 

We also have that dL = 6j — 0^ and so, using (P = 0, it follows that (iOj = 
dQ^- This means that there exists a unique 2-section D,l — — rf0j = —dQ^, that 
will be called the Poincare-Cartan 2-section. This 2-section will be important 
for studying symplecticity of the discrete Euler-Lagrange equations. 

Proposition 4.2. If X and Y are sections of the Lie algebroid AG then 

j7^(x(i.o), yd.")) = , = 0, 

and 

nL{X^^-'°\Y'^'^^^^)^-X{YL) and F^^'")) = y (Xi). 

Proof. A direct computation proves the result. □ 

Remark 4.3. Remark 4.3. Let g be an element of G such that a{g) = x and 
= y- Suppose that U and V are open subsets of , with x ^ U and y G V, 
and that {Xi} and {Yj} are local bases of T{t) on U and V, respectively. Then, 

{Xf '°\y/°'^^} is a local basis of r(7r^) on the open subset ^-^(C/) n (3-\V). 
Moreover, if we denote by {{X'Y^'°\ {Y3)^°^^'>} the dual basis of 
we have that on the open subset a^^{U) fl l3^^{V) 

Ql = -ML){x^Y^'°) , e+ = )(«■!), 

o 

Finally, we obtain some useful expressions of the Poincare-Cartan 1-sections 
using the Lie algebroid isomorphisms introduced in Theorems 13.31 and 13.41 
We recall that the maps 

A{<^P) : A{Gp*Gc,)'^V(3®GVa A{V'^G) ^ V^iAG) 

A{<^'') : A{Gp*Ga)^V(3®GVa A{V"G) ^ V"{AG) 

given by ()3.1ip and ()3.13p are Lie algebroid isomorphisms over the identity of G. 
Moreover, if (we(/3(g)), ^g) e VgiAG) then, from ((OT|) . it follows that 

Ag{^^)-^{v^^p(g)),Zg) = {Z g - (T.^ff {g))l g) (v e(f3 (g))) : iTe{f3{g))lg)iVe{f3{g))))- (4.2) 

On the other hand, if {v^(^a{h))7 Zh) G 'P^{AG) then, using (|3.13p . we deduce that 
Ah{^°')~'^{v^(a(h)),Zh) = (r,(Q(,i))(r^oi)(w,(„(;j))), Z,i-T,(Q(/i))(r,iOi)(w,(„(,j)))). (4.3) 
Now, we introduce the sections 0^ e r((T")*) and 9^ G r((T'^)*) given by 

e2 = (A(<i>")-Md)*(ez), = (A(<i>'^)-Md)*(e+). (4.4) 
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Using (ji?^ and (031), we obtain that 

QLih){Ve(^a(h)),Zh) = ~V^(a(h)){LorhOi), (4.5) 

^L{9){Vc(P(g)),Zg) = V^(p(g)){L O Ig), (4.6) 

for («e(a(/.)),^h) e V^XAG) and («,(;3(g)) , Zg) e T^^^l^G). 

4.2.1. Poincare-Cartan 1-sections: variational motivation. Now, we follow a vari- 
ational procedure to construct the 1-sections 9 J and 6^. We begin by calculating 
the extremals of SL for variations that do not fix the point g £ G. For it, we 
consider the manifold 

= {(51, . . . ,.gAr) e /(.9fc,5fe+i) e G2 for each fc, 1 < A: < iV - 1}. 

If c : (— e, e) ^ is a curve in and c(0) = (gi, . . . , qn) then there exist -I- 1 
curves hk : (-£,e) ^ a-'^{l3{gk)), for < A: < iV, with /ifc(O) = e(/3(gfc)) and 
5o = Sf^i such that 

c{t) = {hQ^{t)gihi{t), h^^{t)g2h2{t), h~\{t)gN-ihN-i{t), h~\{t)gNhN{t)) 

for t G (— e,e). Thus, the tangent space to at (171, . . . , (Jat) may be identified 
with the vector space Api^g^-^G x Ap^g^^G x • • • x Ap(^g^)G, that is, 

^(si,S2,...,9iv)C^ = {(wo,i'i, . . . ,'yAr) / wfe e A^fcG,Xfe = (3{gk),0 <k<N}. 
Now, proceeding as in Section [4.11 we introduce the action sum 

N 

5L:C^^M, ^L(gi,...,gw) = ^L(.gfe). 

Then, 

- 5L(c(t)) = ^ [d°(Lo/,J(e(xfc))(«fe)+d°(ior3,^,oi)(e(a:fc+i))(«fc)] 

+d°{L o o j)(e(a:o))(wo) + d°{L o /<,„)(e(a:w))(«^). (4.7) 

Therefore, if Xq, . . . ,Xjv are sections of r : ylG M satisfying, Xk{xk) — Vk, for 
all k, we have that 

d 



dt 



SL{c{t)) = \Xk{gk){L) - Xk{gk+i){L)\ ~~ xl{gi){L) -f Xn{9n){L) 

" k=l 

N-l 

= E(^DELi(5fc,5^-+i))K) + eZ(5i)(X^'^°^(5i)) + e+(g^)(x(;'^)(5^)). 
fc=i 

Note that it is in the last two terms (that arise from the boundary variations) where 
appear the Poincare-Cartan 1-sections. 

4.3. Discrete Lagrangian evolution operator. We say that a differentiable 
mapping ^ : G — > G is a discrete flow or a discrete Lagrangian evolution 
operator for L if it verifies the following properties: 

- graph(^) C G2, that is, (5,^(5)) £ G2, V.g £ G (^ is a second order 
operator). 

- {g,£,{g)) is a solution of the discrete Eulcr-Lagrange equations, for all g £ 
G, that is, (i^DEL-^) (5=^(5)) = fo^ all 9 e G. 
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In such a case 

d°{Lolg + Lor^^g^^t){e{P{g)))^A,,,,G^O, for all g S G (4.8) 
or, in other terms, 

X{g){L)-X{^{g)){L)^Q (4.9) 

for every section X of AG and every g ^ G. 

Now, we define the prolongation V^S, : VP (Bg Vet — > VjS (Bo Va of the second 
order operator ^ : G — > G as follows: 

A{^°')-^ o{Id,T^)oA{^^) (4.10) 

with A($") and A($^) the isomorphisms defined in Theorems 13.31 and 13.41 and 
{Id,T^) : V'^iAG) V^iAG) the map given by 

{Id,TO{v,Wg)),Xg) = {v,ipig)),{TgO{Xg)), for {v,^p^g)),Xg) e V^{AG). 

Since the pair {{Id, T^), ^) is a Lie algebroid morphism between the Lie algebroids 
V'^{AG) — y G and ^^{AG) — > G then the pair (P^^,0 is also a Lie algebroid 
morphism 

V(3 ®G Va — ^ V(3 ®G Va 



G ^G 




From the definition of V^S,, we deduce that 

ViaX,,Y,) = ((Tg(r,^(,)Oz))(y,),(rgO(^g)+(T;C)(rg)-rg(r,^(,)Oz)(yg)) (4.11) 
for all {Xg,Yg) e Vg(3 ® Vga. Moreover, from ([2?T0l) and (|4TT|) . we obtain that 

V^i{X{g), Y{g)) = {~Y{i{g)), {Tgi){X{g) + Y {g)) + F (^(g))) (4.12) 
for all X, Y sections of AG. 

4.4. Preservation of Poincare-Cartan sections. The following result explains 
the sense in which the discrete Lagrange evolution operator preserves the Poincare- 
Cartan 2-section. 

Theorem 4.4. Let L : G — > R 6e a discrete Lagrangian on a Lie groupoid G. 
Then: 

(i) The map ^ is a discrete Lagrangian evolution operator for L if and only if 

(P-C,0*©Z = 0^ 

(ii) The map ^ is a discrete Lagrangian evolution operator for L if and only if 

{v^uyei-Ql = dL. 

(iii) If is discrete Lagrangian evolution operator then {V^^,S,)*^l — ^l- 
Proof. From ()4.4p . it follows 

(A($"),w)*(e2) = e^, (A(<i>^),M)*(e^) = e+. (4.i3) 

On the other hand, if (^^(^(g)), Xg) G 'P^{AG) then, using (|4.5p and (|4.6p . we have 
that 

{((Id,rO,0*(e2)}(3)(«.(/3(s)),^s) = -^.(P(9)){Lor^^g)0^) 

and 

^l{9){Ve{p{g)),Xg) = We(/3(g))(ioZg). 
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Thus, {{ld,TO,0*Ql = ©L if and only if ^ is a discrete Lagrangian evolution 
operator for L. Therefore, using this fact and (j4.13p . we prove (i). 

The second property follows from (i) by taking into account that dL = O J — 0^ . 
Finally, (iii) follows using (ii) and the fact that {V^S,, £,) is a Lie algebroid morphism. 

□ 

Remark 4.5. Now, we present a proof of the preservation of the Poincare-Cartan 2- 
section using variational arguments. Given a discrete Lagrangian evolution operator 
f : G — > G for i, we may consider the function iS^L : G — > M given by 

iS^L){g)^Lig) + mg)), for .g G G. 

If d is the differential on the Lie algebroid V/3 (Bq Vet G and X, Y are sections 
of r : AG ^ Af then, using ([4A|) . (gj]) and (|4?T2l) . we obtain that 

d{S^L){g){X{g), Y (g)) = X(g)L + Y {g)L + {TS(X{g))L + {T^iY {g))L 

= X{g)L + Y{g)L - Y{C{g))L + {TS(X{g))L 

+ {TSiy{9))L+Y{^{g))L 

= -Ql{g)(X{g),Y{g))+[{V^^,irQV\ {9){X{g),Y{g)). 

This implies that 

Thus, we conclude that = f^L- o 

4.5. Lie groupoid morphisms and reduction. Let ($,$o) be a Lie groupoid 
morphism between the Lie groupoids G =4 M and G' =t M'. The prolongation 
V^^ : Vf3 ®G Va — > Vf3' ©g' Va' of the morphism ($, $o) is defined by 

VmV,W) ^ iTg^V),Tg<P{W)) (4.14) 

for every (V, W) € Vg/3 © VgO;. It is easy to see that {V^^, $) is a morphism of Lie 
algebroids. 

Theorem 4.6. Let ($,<f>o) be a morphism of Lie groupoids from G ^ M to G' ^ 
M'. Let L and L' be discrete Lagrangian functions on G and G' , respectively, 
related by L = L' o ^ . Then: 

(i) for every (g, h) e G2 and every v G A^(-g-|G we have that 

DuELL{g,h){v) = DuELL'mg),^h)){Afj^g~,^v)). (4.15) 

(ii) (P-$,ci>)*e+ =6+ 

(iii) {v^<f,^)*ei, = ei, 

(iv) {p-<^,<i>)*nL' = ^L. 

Proof. To prove the first we notice that, if ($, $0) is a morphism of Lie groupoids, 
then we have that ^ o Ig = o $ and $ o r/i = r$(-;j-) o $, from where we get 

DoEhLig, h){v) = Tlg{v)L + Trh{Ti{v))L 

= Tlg{v){L' o $) + Trh{Ti{v)){L' o $) 

= T<^{Tlg(v))L' + T<^{Trh{Ti(v)))L' 

= Th(g){T^{v))L' + Tr^^h){T^Ttiv)))L' 

= T/$(<,)(T$(«))L' + Tr$(^)(Tz'(T$(t-)))L' 

= DuELL'mg),^h)){Ap^g)^v)), 

where we have also used that z' o <[) = $ o i and A^(g)<i>(ti) = T$(ti). 
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For the proof of the second, we have that 

{{V^^,^ret,ig),iV,W)) - {Q+mg)),iT,^V),Tg^W))) 

= {{Tg<S>)iW))L' ^WL= (e+(g), {V, W)), 

for every (V, W) £ T^gG. The proof of the third is similar to the second, and finaUy, 
for the proof of (iv) we just take the differential in (ii). □ 

As an immediate consequence of the above theorem we have that 

Corollary 4.7. Let ($,(l>o) be a morphism of Lie groupoids from G ^ M to 
G' =1 M' and suppose that (g, h) G 6*2. 

(i) // ($((7), is a solution of the discrete Euler- Lagrange equations for 
L' — Lo<i), then (g, h) is a solution of the discrete Euler- Lagrange equations 
for L. 

(ii) //$ is a submersion then (5, h) is a solution of the discrete Euler- Lagrange 
equations for L if and only if (^{g),^{h)) is a solution of the discrete 
Euler- Lagrange equations for L' . 

(iii) // $ is an immersion, then (g, h) is a solution of the discrete Euler- 
Lagrange equations for L if and only if DY)EhL{^{g),^{h)) vanishes over 
Im(A/3(g)<i>). 

The case when $ is an inmersion may be useful to modelize holonomic mechanics 
on Lie groupoids, which is an imprescindible tool for explicitely construct geometric 
integrators (see [51 \7T\). 

The particular case when $ is a submersion is relevant for reduction (see Section 
15.51 in this paper) . 

4.6. Discrete Legendre transformations. Given a Lagrangian L : G — > K we 
define, just as the standard case [21], two discrete Legendre transformations 

¥-L : G — > A*G and ¥+L : G — > A*G as follows 

L){h){v^(^a{h))) = -Ve{a(h)){LorhOi), ior V^(a{h)) & Aa{h)G, 

(F+L)(.g)(t;,(^(g))) = ^.(^(g)) (i O /g), for G A^(g)G. 

Remark 4.8. Note that {¥~L){h) S ^a(h)G- Furthermore, if U is an open subset 
of M such that a{h) G U and {Xi} is a local basis of T{t) on U then 

¥-L = xl{L){X'' oa), 

on a^^{U), where {X^} is the dual basis of {X^}. In a similar way, if V is an open 
subset of M such that P{g) G V and {Yj} is a local basis of T{t) on V then 

¥+L = Yj{L)iY^ op), 

onp-'^{V). o 

Next, we consider the prolongation : (AG) A*G oi the Lie algebroid 
T : AG ^ M over the fibration r* : A*G M, that is, 

VlliAG) = {(«,.(,.), X,.) eA,.(..)GxT,.(A*G)/(T,.(,.)/3)(i;,.(,.)) 
= (T,.T*)(X„0} 

for V* G A*G. Then, we may introduce the canonical section 8 of the vector bundle 
(r^*)* : {V' AG)* A*G as follows: 

&{v*){Vr.i..),X,.) = !'*(«,.(„.)), (4.16) 

for V* G A*G and (wt-.(^.), X„.) G V^-, {AG). O is called the Liouville section. 
Moreover, we define the canonical symplectic section il associated with AG by 
fl = -dQ, where d is the differential on the Lie algebroid r'^* : V^* {AG) A*G. 
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It is easy to prove that is nondegenerate and closed, that is, it is a symplectic 
section of P^' (AG) (see [HI [23]). 

Now, let V^F^L be the prolongation of W^L defined by 

V''¥-L={ld,T¥-L)oA{<S>°') ■.r^G = V/3®GVa — ^V'iAG), (4.17) 

where A($") : VG = F/3 ®g Va V^iAG) is the Lie algebroid isomorphism 
(over the identity of G) defined by (|XT3)l and {Id,T¥-L) : V"{AG) {AG) is 

the map given by 

for {v^(a{h)),Xh) e V^{AG). Since the pair {{Id,T¥^ L),W^ L) is a morphism be- 
tween the Lie algebroids ViAG) G and (AG) A*G, we deduce that 
{V^¥~L,¥~L) is also a morphism between the Lie algebroids V^G = VP(BGVct 
G and i^G) — > A*G. The following diagram illustrates the above situation: 

V/3 ©G Va '^^^ ^ V^' (AG) 




M 



The prolongation V^¥ L can be explicitly written as 
Vl¥'L{Xh,Yh) = inilorh-^){Xh),iTh¥-L)iX^,) + {Th¥- L){Yh)) , (4.18) 
ioiheG and {Xh, Yh) G © V^a. 

Proposition 4.9. If Q is the Liouville section of the vector bundle {V^ {AG))* 
A*G and Vl — ~dO is the canonical symplectic section of A^lV^ (AG))* A*G 
then 

{P^(¥-L),¥-L)*e = e^, {r^i¥-L),¥-L)*n = Ql- 

Proof Let 92 be the section of (r")* : V^iAG)* G defined by Then, 
from (|4.5p and (|4.16p . we deduce that 

{{Id,T¥-L),¥-L)*e = Q2. 

Thus, using (|4.4p . we obtain that 

{P''¥-L,¥-L)*Q = 6^. 

Therefore, since the pair {P'^¥^L, ¥^L) is a Lie algebroid morphism, it follows that 

{r''¥-L,¥-L)*n = VLl- 

□ 

Now, we consider the prolongation V'^¥^L of F+L defined by 

^¥+1 = {Id, T¥+L) o A{¥^) : V^G = V/S ®g Va — > {AG), (4.19) 

where A{^l^) : VG = VP ©g Va 'P^{AG) is the Lie algebroid isomorphism 
(over the identity of G) defined by (IXTTI) and {Id, T¥+L) : V'^{AG) V' {AG) is 
the map given by 

{Id,T¥+L){v,(p(,)),Xg) = {v,^p(,)),{Tg¥+L){Xg)), 
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for (we(/3(g)), Xg) e P^{AG). As above, the pair (P'^F+L, F+i) is a morphism be- 
tween the Lie algebroids V^G = VP QcVa ^ G and P''' (AG) A*G and the 
following diagram illustrates the situation 

^^""^ - -'{AG) 



^A*G 



We also have: 

Proposition 4.10. 7/0 is the Liouville section of the vector bundle {V^ {AG))* — > 
A*G and O = —dQ is the canonical symplectic section of A^{V'^* {AG))* —>■ A*G 
then 

(7'^(F+L),F+L)*e = et, {P^¥+L),¥+L)*n ^ Ql- 

Remark 4.11. (?) If ^ : G — > G is a smooth map then ^ is a discrete Lagrangian 
evolution operator for L if and only if F~L o ^ = F+L. 

{ii) If {g, h) e G2 we have that 

(^DELi)(5, h) = ¥+L{g) - ¥-L{h). (4.20) 



4.7. Discrete regular Lagrangians. First of all, we will introduce the notion of 

a discrete regular Lagrangian. 

Definition 4.12. A Lagrangian L : G on a Lie groupoid G is said to be regular 
if the Poincare-Cartan 2-section SIl is symplectic on the Lie algebroid P'^G = 
Vl3 ®G Va G. 

Next, we will obtain necessary and sufficient conditions for a discrete Lagrangian 

on a Lie groupoid to be regular. 

Theorem 4.13. Let L : G ^M. be a Lagrangian function. Then: 
a) The following conditions are equivalent: 

(i) L is regular. 

(ii) The Legendre transformation ¥~L is a local diffeomorphism. 

(iii) The Legendre transformation F+L is a local diffeomorphism. 

h) If L : G ^ W is regular and {go,ho) G G2 is a solution of the discrete 
Euler- Lagrange equations for L then there exist two open subsets Uq and Vq of G, 
with go e Uq and ho G Vo, and there exists a (local) discrete Lagrangian evolution 
operator £,l ■ Uq ^ Vo such that: 

(i) Cl(.9o) = ho, 

(ii) is a diffeomorphism and 

(iii) is unique, that is, if Uq is an open subset of G, with go e Uq and : 
Uq ^ G is a (local) discrete Lagrangian evolution operator then Cz,|i7on(7^ ~ 
?L|c/onc/^- 



Proof 



a) 



First we will deduce the equivalence of the three conditions 



(i)^(ii) If h e G, we need to prove that Th{¥-L) : ThG — > Tr-L(h)A*G is a 
linear isomorphism. Assume that there exists g T^G such that Th(¥^ L){Yh) = 
0. Since r* o F'L = a, then {Tha){Yh) = 0, that is, Yu e Vna. 
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Therefore, (0,,, Yh) e Vhl3®Vha and, from (liT^ . wc have that Vl{¥- L){Oh, Y,,) = 
0. Moreover, lvi{V- L)YVt{¥- L{h)) = ^^^(/i) and f7(F-L(/i)) and are non- 

degenerate. Therefore, we deduce that VJ^iW^ L) is a hnear isomorphism. This 
imphes that Yu = 0. This proves that ThiV' L) : TuG Tf-L(h){A*G) is a hnear 



isomorphism. In the same way we deduce (i)=>(iii) 



(ii)=>(i) Wc will assume that F+L is a local diffconiorphism, so that 



Vl^+L : VIG ^ V,(3 © V,a V^l^^^^iAG) 

is a linear isomorphism, for all 5 G G. 

On the other hand, if is the canonical symplectic section of the vector bundle 
h^{V^'{AG))* A*G then, from Proposition QUI we deduce that 

(p;F+L)*(f7(F+i(g))) = ^L{g). 

Thus, since Vl{¥^ L{g)) is nondegenerate, we conclude that 0^(17) is also nondegen- 



erate, for all g £ G. Using the same arguments we deduce (iii)=>(i) 



b) Using Remark 4.11, we have that 

(F+i)(go) = (¥-L)(ho) = Mo e A*G. 

Thus, from the first part of this theorem, it follows that there exit two open subsets 
Uo and Vq of G, with go £ Un and ho G Vo, and an open subset Wo of A*G such 
that Ho & Wo and 

¥+L : Uo Wo, ¥-L : Vo Wo 
are diffeomorphisms. Therefore, using Remark 4.11, we deduce that 

a = [(F-L)-i o (F+L)]|^„ -.Uo^Vo 

is a (local) discrete Lagrangian evolution operator. Moreover, it is clear that 
Ci(5o) = ho and, from the first part of this theorem, we have that S,l is a dif- 
feomorphism. 

Finally, ii Uq is an open subset of G, with go Uo, and : ?7q ^ G is another 
(local) discrete Lagrangian evolution operator then ^^jj^p^iji^ : UoCiUq ^ G is also 
a (local) discrete Lagrangian evolution operator. Consequently, using Remark 4.11, 
we conclude that 

□ 

Remark 4.14. Using Remark 4.3, we deduce that the Lagrangian L is regular 
if and only if for every g € G and every local basis {Xi} (respectively, {Yj}) of 
T{t) on an open subset U (respectively, V) of M such that a{g) G U (respectively, 
Pig) G V) we have that the matrix Xi{Yj{L)) is regular on a~'^{U) n (3~^{V). o 

Let L : G ^ M be a regular discrete Lagrangian on G. If / : G ^ M is a real 
G°°-function on G then, using Theorem I4.13( it follows that there exists a unique 
^/ G r(7r^) such that 

lij^L^df, (4.21) 

d being the differential of the Lie algebroid tt^ : V^G = V (3 ®g Vot G. ^/ is 
called the Hamiltonian section associated to f with respect to Q.^. 

Now, one may introduce a bracket of real functions on G as follows: 

{■,-}l:C°°{G)xC°^{G)^C°^{G), {f,g}L = -^L{^f,Q. (4.22) 
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Note that, from (|4.2ip and Propositions 14.91 and 14.101 we obtain that 

(P-F±L) o e/oF±L = ° F±L, (4.23) 

for / e C°°{A*G), where V^F^L : VG = VP ©g Va V' (AG) is the pro- 
longation of F^L (see Section |4!6|) and ^/ is the Hamiltonian section associated to 
the real function / on A*G with respect to the canonical symplectic section D, on 
A^iV' (AG))* A*G, that is, i^^M = df. 

On the other hand, we consider the canonical linear Poisson bracket {•,•} : 
G'^{A*G) X C°°(A*G) ^ G°°{A*G) on A*G defined by (see [HI) 

{/,g} = -f^(C/,C,), for/,geC°°(A*G). (4.24) 

We have that (see [H]) 

Moreover, from and Propositions li^ and BTTUl we deduce that 

{/ o F±L, g o F±L}i = {/, g} o F±L. 
Using the above facts, we may prove the following result. 

Proposition 4.15. Let L : G ^ M be a regular discrete Lagrangian. 

(i) The Hamiltonian sections with respect to fl^ form a Lie suhalgebra of the 

Lie algebra (T{T:^),l-,-f^^)- 
ill) The Lie groupoid G endowed with the bracket {•, -II is a Poisson manifold, 
that is, {-, -ji is skew- symmetric, it is a derivation in each argument with 
respect to the usual product of functions and it satisfies the Jacobi identity. 
(iii) The Legendre transformations F^L : G ^ A*G are local Poisson isomor- 
phisms. 

4.8. Discrete Hamiltonian evolution operator. Let i : G ^ R be a regular 
Lagrangian and assume, without the loss of generality, that the Legendre transfor- 
mations F+L and F^ L are global diffeomorphisms. Then, — (F^i)^^ o (F+L) 
is the discrete Euler-Lagrange evolution operator and one may define the discrete 
Hamiltonian evolution operator, : A*G ^ A*G, by 

a =F+Lo^io(F+i)-i . (4.25) 

From Remark 14.111 we have the following alternative definitions 

a - F-L o o (F-i)-i, - F+L o (F-L)^i 

of the discrete Hamiltonian evolution operator. The following commutative diagram 
illustrates the situation 




A*G 



A*G 



A*G 



Define the prolongation : (AG) V^' {AG) of |l by 

V'''Il = V^F+L o V^S,L o {V^F+L)-\ 
or, ahernatively (see ((iTII)) . (|iT7)) and dHH)), 

P^'iL = V^F+L o {P^F-L)-\ V'^L = VF-L o V^(,l o {VF-L^^. (4.26) 
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Proposition 4.16. //8 is the Liouville section of the vector bundle [V^ {AG))* 
A*G and Vl — —dO is the canonical symplectic section of A^{V^ (AG))* — > A*G 
then 

[v^'iL-.^Lre = e + d(Lo (F-L)-i), = n. 

Moreover, is a Poisson morphism for the canonical Poisson bracket on A*G. 

Proof. The result follows using (|4.25|) . (|4.26p and Theorem 14.41 and Propositions 
lOand liTSl □ 

4.9. Noether's theorem. Recall that classical Noether's theorem states that a 
continuous symmetry of a Lagrangian leads to constants of the motion. In this 
section, we prove a discrete version of Noether's theorem, i.e., a theorem relating 
invariance of the discrete Lagrangian under some transformation with the existence 
of constants of the motion. 

Definition 4.17. A section X of AG is said to be a Noether's symmetry of the 
Lagrangian L if there exists a function f G C°° {M) such that 

dL{X^^'^'>)^p*f-a*f. 

In this case, L is said to be quasi-invariant under X. 

When = -XL + XL^O, we will say that L is invariant under X 

or that X is an infinitesimal symmetry of the discrete Lagrangian L. 

Remark 4.18. The infinitesimal invariance of the Lagrangian corresponds to a 
finite invariance property as follows. Let $s the flow of X and 7(5) = ^s{^{x)) 
be its integral curve with 7(0) = e{x), where x = P{g). Then, the integral curve 
of X at (7 is s t-^ r^(s)g = gi{s), and the integral curve of —X through e(x) is 
s ^ 7(5)""'^. On the other hand, if [h, h') G G2 and Yh G VhP, Z^t e Vh'OL then 

Using the above facts, we deduce that the integral curve \i of the vector field ~X^X 
on G satisfying /Lt(0) = g is 

= 7(5)^^97(5)7 for aU s. 
Thus, the invariance of the Lagrangian may be written as 
i(7(s)-'57(5)) = L{3), for all s. 

o 

If i : G ^ R is a regular discrete Lagrangian, by a constant of the motion we 
mean a function F invariant under the discrete Euler-Lagrange evolution operator 
a, that is, f =F. 

Theorem 4.19 (Discrete Noether's theorem). If X is a Noether symmetry of a 
discrete Lagrangian L, then the function F = Q~j^{X'^^'^^) — a* f is a constant of 
the motion for the discrete dynamics defined by L. 

Proof Wc first notice that e^(X(i'i)) = XL so that the function is = XL - 
a*f. 

If the Lagrangian L is quasi-invariant under X and g is a point in G, then 
-Xig)iL) + X{g){L) - /(/3(g)) - /(a(g)), 

so that 

X{g){L) = X{g){L) + f{P{g)) - fHg))- 
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We substrate X {^L(g)){L) to both sides of the above expression, so that 

Xig){L) X{U9)m = [X{g){L) f{a{g))] [X{i{g)){L) f{a{^L{9))] 
^F{g)-F{iL{9)), 

from where the resuh immediately foUows using (j4.9p . □ 
Proposition 4.20. If X is a Noether symmetry of the discrete Lagrangian L then 

£^(1,1)9^ =d(a7). (4.27) 

Thus, if L is regular, the complete lift X''^'^^ is a Hamiltonian section with 
Hamiltonian function F = — a* f , i.e. z^(i,i)f^L = dF . 

Proof Indeed, if dL(X(i'i)) = (3* f - a* f and F is a section of AG, we have that 
(see Proposition imi), 

(/:xa..ez)(i^(^'°^) = -f^i(x(i'i),y(i'0)) + d(z;,a.,ez)(r(i'°)) 

= -Y{XL) + Y(XL)^Y{a*f-P*f) 
= d{a* f){Y^'''^y). 
On the other hand, using (|2.9p and Proposition l4.21 we deduce that 

(/:xa..ez)(r(°-^^) = -f7L(x(i-i),r(°-i)) + d(^^aa,e^)(r("-i)) 

= -X{YL) + Y(XL) = \Y, X]{L) = = d{a* f){Y^'''^^). 
Thus, KT!\ holds. From (14271) . it follows that 

which completes the proof. □ 
We also have 

Proposition 4.21. The vector space of Noether symmetries of the Lagrangian 
L : G —I- R is a Lie subalgebra of Lie algebra {T{t), |., •]). 

Proof. Suppose that X and Y are Noether symmetries of L and that 

dL(X(i'i)) = -XL + XL = f3*f - a*f, (4.28) 

di(y(i.i)) = - YL +YL = f3*g - a*g, (4.29) 
with /,5 e C°°(M). Then, using ([211), ((321) and jHSl), we have that 

dL{lX, rl^i'i)) = X{YL) - y(XL) + X{YL) - Y{XL). (4.30) 
On the other hand, from (j^ . (P?^ and we deduce that 

X(yL) = X(ri) - a*(p(X)(.g)), y (XL) = YiXL) - a*ipiY){f)), 
X{YL) = X{YL) + [3*{p{X){g)), Y{XL) = y(XL) + /3* {p{Y){f)). 

Thus, using ((T^ and (IQU]) . we obtain that 

dL(lX, rjd'i)) ^ f3*{p{X){g) ~ p{Y){f)) - a*{p{X){g) p{Y){f )). 

Therefore, {X, Yj is a Noether symmetry of L. □ 

Remark 4.22. If L : G R is a regular discrete Lagrangian then, from Proposi- 
tions 14.201 and 14.211 it follows that the complete lifts of Noether symmetries of L 
are a Lie subalgebra of the Lie algebra of Hamiltonian sections with respect to Hl. 
o 
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5. Examples 

5.1. Pair or Banal groupoid. We consider the pair (banal) groupoid G — MxM, 
where the structural maps are 

a{x,y) = x, (3{x,y)^y, e{x) ^ {x,x), i{x,y) ^ {y,x), 
m{{x,y), {y,z)) = {x,z). 

We know that the Lie algebroid of G is isomorphic to the standard Lie algebroid 
TM ■■ TM M and the map 

^ ■.AG = V,i^M)a ^ TM, (0„ w,) e T,M x T^M ^ «',(0„ v,) =v,, for a; e M, 
induces an isomorphism (over the identity of M) between AG and TM. If X is a 
section of tm ■ AG ~ TM M, that is, X is a vector field on M then X and X 
are the vector fields on M x M given by 

Xix, y) - {~X{x),%) e T,M X TyM and X{x, y) = (0,, G T,M x TyM, 

for (x, y) ^ M X M. On the other hand, if (x, y) e M x M we have that the map 
Vl^'y^G = l/(,,j,)/3 ® V(^,y)a ^ r(,,y)(M X M) ^ T,M x TyM, 

induces an isomorphism (over the identity of M x M) between the Lie algebroids 
T^TM . -prMQ ^vpcSaVa^ G ^ M x M and T(mxa/) : T{M x M) ^ M x M. 

Now, given a discrete Lagrangian L : M x M ^ R then the discrete Euler- 
lagrange equations for L are: 

X{x,y){L)-X{y,z){L) = 0, for all X e X(M), (5.1) 

which are equivalent to the classical discrete Euler-Lagrange equations 

D2L{x,y) + DiL{y,z)^0 

(see, for instance, [21] )■ The Poincare-Cartan 1-sections 8^ and 0j on n'^'^' : 
■prMQ ^ X M) ^G = M X M are the 1-forms on M x M defined by 

'dLi^^y)i'"x,Vv) = -Vx{L), Ql{x,y){va;,Vy) = Vy{L), 

for {x,y) e M x M and ivx,Vy) G T^M x T^^Af ~ T(^,y){M x M). 

In addition, if ^ : G = M x M ^ G = M x M is a discrete Lagrangian evolution 
operator then the prolongation of ^ 

-pTM^ . -prMQ ^ X M) -P^"G ~ T(M X M) 

is just the tangent map to ^ and, thus, we have that 

= —dQ^ = — rfOj being the Poincare-Cartan 2-form on M x M. The Legendre 
transformations ¥- L : G = M x M ^ A*G :^ T*M and¥+L : G = M x M 
A*G ^ T*M associated with L are the maps given by 

F^L(x, y) = -DiL{x, y) e T^M, ¥+L{x, y) = D2L{x, y) e T;M 

for {x, y) G Af X M. The Lagrangian L is regular if and only if the matrix 

is regular. Finally, a Noether symmetry is a vector field X on Af such that 

D,Lix,y){X{x)) + D^Lix, y)iX{y)) = f{y) - f{x), 

for {x,y) e M X M, where / : Af ^ M is a real G°°-function on M. If X is a 
Noether symmetry then 

X F{x) = D.Lix, y)iXix)) - f{x) 

is a constant of the motion. 
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In conclusion, we recover all the geometrical formulation of the classical discrete 
Mechanics on the discrete state space M x AI (see, for instance, 21j). 

5.2. Lie groups. We consider a Lie group G as a groupoid over one point M = {e}, 
the identity element of G. The structural maps are 

= e, f3{g)^t, e(e) = e, i{g)=g~'^, m{g,h)=gh, iorg,heG. 

The Lie algebroid associated with G is just the Lie algebra g — T^G of G. Given 
^ e we have the left and right invariant vector fields: 

a9) = (Tcg(e), ti9)^iT,r,m, for^eC. 

Thus, given a Lagrangian L : G — > M its discrete Euler-Lagrange equations are: 

(TJg^miL) - (T.rg^^^miL) = 0, for all C G fl and 5fe,5fc+i e G, 

or, (r*^dL){c) = (r*^^^dL)(e). Denote by = (^gfc^-^)(^) then the discrete Euler- 
Lagrange equations are written as 

Hk+i ^ Ad*g^fik, (5.2) 

where Ad : Gx g — > g is the adjoint action of G on g. These equations are known 
as the discrete Lie-Poisson equations (see [T| fl ^ 1 ^ ). 

Finally, an infinitesimal symmetry of L is an element ^ S fl such that {Tclg){(,) (L) = 
(re?'g)(^)(L), and then the associated constant of the motion is F((7) = {TJg){£^){L) — 
(Terg)(^)(i). Observe that all the Noether's symmetries are infinitesimal symme- 
tries of L. 

5.3. Transformation or action Lie groupoid. Let H he a, Lie group and ■ 
M X H ^ M, {x, h) <E M X H ^ xh, a right action of H on M. As we know, H 
is a Lie groupoid over the identity element e oi H and we will denote by a, (3, e, m 
and i the structural maps of H . If tt : M — > {e} is the constant map then is clear 
that the space 

x^H ^ {(.T, h) eM X H/tt{x) = a{h)} 

is the cartesian product G — M x H and that ■ : M x H ^ M induces an action 
of the Lie groupoid H over the map tt : M ^ {e} in the sense of Section 12.21 
(see Example 6 in Section [^72]) . Thus, we may consider the action Lie groupoid 
G — M X H over M with structural maps given by 

a^{x,h)^x, (3^{x,h)^xh, e^(a;) = (x, e), , s 

m7r((a;, h), {xh, h')) = {x, hh'), h) — {xh, h^^). 

Now, let \) — T^H be the Lie algebra of H and <I> : f) X(M) the map given by 

^{v) = »/M, for ?/ e [), 

where rjM is the infinitesimal generator of the action ■ : M x H ^ M corresponding 
to rj. Then, $ defines an action of the Lie algebroid (} — > {a point} over the 
projection tt : M — > {a point} and the corresponding action Lie algebroid pri : 
M X [) — > M is just the Lie algebroid of G = M x H (see Example 6 in Section [2?2l) . 

We have that T{pri) = {fj : M t)/^ is smooth } and that the Lie algebroid 
structure (|-, •]$, p$) on pri : M x H M is given by 

[fl,ill^{x) = [f]{x),il{x)] + {fi{x))M{x){il)-{fi{x))M{x){fi), p<s>{'n){x) = {fi{x))M{x), 

for Tj,fi€ r(pri) and x G M. Here, [•, •] denotes the Lie bracket of f). 

If {x, h) G G = M X H then the left-translation l(x.h) ■ ct^^{xh) a^^(a:;) and 
the right-translation r^^j^'f : l3~^{x) l3~^{xh) are given 

li.,h){xh, h') = {x, hh'), r^,j,^{x{h')-\h') = {x{h')-\h'h). (5.4) 
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Now, if 77 G f) then 77 defines a constant section : M — > () of pri : M xi) ^ M 
and, using (j2.6p . ()2.7p . (15. 3p and (|5.4p . we have that the left-invariant and the 
right-invariant vector fields and C,,, respectively, on M x H are defined by 

Cjj{x,h) = {-riM{x),Ti'{h)), C^(a;, ft.) = (Oo;, (5-5) 

for (a;, h) e G = M x H. 

Note that if {rji} is a basis of [} then {C^^} is a global basis of r(pri). 

Next, suppose that L : G = M x H ^ is a. Lagrangian function and for every 
h ^ H (resp., x G M) we will denote by Lh (resp., L^) the real function on M (resp., 
on H) given by Lh{y) = L{y^ h) (resp., Lx{h') = L{x, h')). Then, a composable pair 
((x, /ife), (xhk, hk+i)) G G2 is a solution of the discrete Euler-Lagrange equations 
for L ii ^ _^ 

Cn{x,hk){L) - C r,{xhk,hk+i){L) = 0, for aU 77 € [), 
or, in other terms (see (|5.5p ) 

{{T.lhM}{L,) - {{T^rn,^,)m{L.h,) + m{xhk){Lh,^,) = 0, for aU ,7 e t). 

As in the case of Lie groups, denote by ^kix,hk) — d[Lx o rh^^){e)- Then, the 
discrete Euler-Lagrange equations for L arc written as 

^j.k+i{xhk,hk+i) = Adl^fikix, hk) + d{Lh^^^ o {{xhk)-)){e), 

where (xhk)- : H —t M is the map defined by 

{xhk) ■ (h) = x{hkh), for h € H. 

In the particular case when M is the orbit of a S ^ under a representation of G 
on a real vector space V, the resultant equations were obtained by Bobenko and 
Suris, see [HE], and they were called the discrete Euler-Poincare equations. 
Finally, an element ^ G t) is an infinitesimal symmetry of L if 

CM{x)iLh) ~ ^{h){L^) + T{h){L^) - f{xh) - fix) 

where / : M — > M is a real C°°-function on M. The associated constant of the 
motion is 

F{x, h) = -Ui{x){Lh) + l{h){L,) - fix), 

for ix, h) gM X H. 

The heavy top. As a concrete example of a system on a transformation Lie groupoid 
we consider a discretization of the heavy top. In the continuous theory [52], the 
configuration manifold is the transformation Lie algebroid t : x so (3) — > 5*^ with 
Lagrangian 

LciT, n)^^n-m~ mglV ■ e, 

where G R'^ ~ so (3) is the angular velocity, F is the direction opposite to the 
gravity and e is a unit vector in the direction from the fixed point to the center of 
mass, all them expressed in a frame fixed to the body. The constants m, g and I 
are respectively the mass of the body, the strength of the gravitational acceleration 
and the distance from the fixed point to the center of mass. The matrix / is the 
inertia tensor of the body. In order to discretize this Lagrangian it is better to 
express it in terms of the matrices f2 G so(3) such that ilw = x w. Then 

ic(r, ^)^i^ T^riVlltf) - mglT ■ e. 

where I ^ \ Tr(/)/3— /. We can define a discrete Lagrangian L : G = S^xSOiS) 
M for the heavy top by 

LiTk,Wk) = Tr(Wfe) - hmglTk ■ e. 
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which is obtained by the rule 0, = R^R « j^Rj^{Rk+i — Rk) = jii^k — h), where 

Wk = RlRk+i- 

The value of the action on an admissible variation is 

= [Tr(Wfce*^) + mglh^Tk ■ e + Tr(Je-*^W^fc+i) + mglh^ie^'^'Tk+i) ■ e] , 

where Tk+i — W'^Tk (since the above pairs must be composable) and K G 5o(3) is 
arbitrary. Taking the derivative at t = and after some straightforward manipula- 
tions we get the DEL equations 

Mk+i - W^MkWk + mglh\r^^ e) = 

where M = WI — IW^ . Finally, in terms of the axial vector 11 in K'^ defined by 
n = M, we can write the equations in the form 

Uk+i = W^FUk + mglh^Tk+i x e. 

Remark 5.1. The above equations are to be solved as follows. From F^, Wfc we 
obtain Ffe+i = VFfeFfe and 11^ from 11^ = Wkl — IW^ . The DEL equation gives 
li-k+i in terms of the above data. Finally we get Wk+i as the solution of the 
equation tik+i = Wk+il — IW^_^^, as in 26J. o 

In the continuous theory, the section X(r) — (r,F) of S"^ x so(3) S*^ is a 
symmetry of the Lagrangian (see [22] )■ We will show next that such a section is 
also a symmetry of the discrete Lagrangian. Indeed, it is easy to see that the left 
and right vector fields associated to X coincide X — X and are both equal to 

X(F, W) = ((F, 0), {W, f VF))) e TG = TS^ x r50(3). 

Thus "-^(X*^^'^)) = so that X is a symmetry of the Lagrangian. In fact it is a 
symmetry of any discrete Lagrangian defined on G = S*^ x 50(3). The associated 
constant of motion is 

(XL){W, F) = Tr(jf VF) = ^ Tr[(W^J - JIF^)f] = -H • F, 

i.e. (minus) the angular momentum in the direction of the vector F. 

5.4. Atiyah or gauge groupoids. Let p : Q — > A/ be a principal G-bundle. A 
discrete connection on p : Q M is a map Ad : Q y. Q ^ G such that 

Ad{gq,hq') = hAd{q,q')g^^ and Adiq,q)^e (5.6) 

for g,h e G and q,q' e Q, t being the identity in the group G (see [TH [H]). 
We remark that a discrete principal connection may be considered as the discrete 
version of an standard (continuous) connection on p : Q — > M . In fact, if Ad '■ 
Q X Q G is such a connection then it induces, in a natural way, a continuous 
connection Ac ■ TQ g defined by 

Ac{Vq) = mq^q)Ad){Oq,Vq), 

for Vq G TqQ. Moreover, if we choose a local trivialization of the principal bundle 
p : Q ^ M to he G X U , where U is an open subset of M then, from (|5.6p . it follows 
that there exists a map A : U x U ^ G such that 

-^diig,x),{g',y)) ^ g'A{x,y)g~\ and A{x,x) = t, 

for {g,x), {g' ,x') € G xU (for more details, see [T^ [T^ l. 

On the other hand, using the discrete connection Ad, one may identify the open 
subset {p-\U) X p-\U))/G ~ ((G x U) x {G x U))/G of the Atiyah groupoid 
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{Q X Q)/G with the product manifold {U xU) x G. Indeed, it is easy to prove that 
the map 

((G X [/) X (G X U))/G -^(UxU)xG, 
[((<?, x), {g', y))] {{x, y),Adiie, x), {g-^g', y))) = {{x, y),g-^g'A{x, yj), 

is bijective. Thus, the restriction to ((G x U) x {G x U))/G of the Lie groupoid 
structure on {Q x Q)/G induces a Lie groupoid structure in ([/ x U) x G with 
source, target and identity section given by 

a:{UxU)xG-^U] {{x,y),g)^x, 
(3:{UxU)xG^U; {{x,y),g)^y, 
e:U^{UxU)xG; x ^ {{x,x),e), 

and with multiphcation m : {{U x U) x G)2 {U x U) x G and inversion i : 
{U xU) X G ^ {U X U) X G defined by 

m{{{x,y),g),{{y,z),h)) = {{x, z), gA{x,yy^hA{y, zy^A{x, z)), 
i{{x,v),g) = {{y,x),A{x,y)g-^A{y,x)). 

The fibre over the point a: e [/ of the Lie algebroid A{{U xU)xG) may be identified 
with the vector space TxU x g. Thus, a section of A{{U x [/) x G) is a pair {X, ^), 
where X is a vector field on U and <^ is a map from U on g. Note that the space 
T{A{{U X U) X G)) is generated by sections of the form (X, 0) and (0, Gj), with 
X £ X{U), S, E g and Gj : J7 — > g being the constant map C^{x) — ^, for all x G U. 
Moreover, an straightforward computation, using ()5.7|) . proves that the vector fields 

iX,0), {X,oi, (0,G^) and (oTQ) on {U x U) x G are given by 

lxj)iix, y),g) = iO,,X{y),{TA(.,y)lgA(..y)-^{{TyA,){X{y)))+ 

-{AdA^,,^yy^{TyAy){X{y))y{g))), 

iX;0){ix,y),g) = (-X(x),0„-(T^(,,,)l,^(,,,)-i((r,A,)(X(a;))))+ 

-{AdAi.,y)-^{TM{X{x))y{g))), 

iO,C^)iix,y),g) = {0,,0y,{AdAi^,y)-iO\9)), 

(0,Gej((a;,y),.g) = (0„0„r(5)), 

for ((x, 2;), g) E {U xU) X G, where Ih : G ^ G denotes the left-translation in G by 
/igG, Ad :Gx0^g is the adjoint action of the Lie group G on g, (respectively, 
rf) is the left-invariant (respectively, right-invariant) vector field on G such that 
77' (e) = T] (respectively, if{t) = i]) and A^ : U ^ G and Ay : U ^ G are the maps 
defined by 

AAy) = Ay{x) ^ A{x,y). 

Now, suppose that L : (Q x (5)/G — > K is a Lagrangian function on the Atiyah 
groupoid {Q x Q)/G. Then, the discrete Euler-Lagrange equations for L are 

lxJ){ {x,y),gk){L) - j X^O) [j y, z), gk+i){L) - 0, 
(0,COiix,y),gk){L) - {Q,C^){{y, z), gk+i){L) = 0, 

with X e X{U), e e and {{{x, y), gu), {{y, z), gk+i)) e {{U x U) x G)2. 
From ()5.8p . it follows that the above equations may be written as 

D2L{{x,y),gk) + DiL((y, z), gk+i) + dfAL[x,y,gk]{y) + 

+ dfAL[y, z, gk+i]{y) + dfX^j[x,y, gk]{y) + dfl^j[y,z,gk + l]{y) = 0, (5.9) 

d{L(x,y,) °lg, ° lAix,y)-^)ii) - C?(i(y,^ J °?'gfc+i)(0 =0, (5.10) 
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where Ig : G G denotes the interior automorphism inG oi g £ G, L(^^ y^ j : G — > R 
is the function L(s^.y^ ){g) = L{x,y,g), and y, g], f\Lj[x,y,g] and fALii^^V^a] 
are the real functions on U given by 

fAL[x,y,g]{y) = L{x,y,gA(x,yy^A{x,y)), 

fALi[x,y,9\{y) = Hx,y,gA{x,y)~^A{y,y)A{x,y)), 

fALi[^,y^9]{y) = L{x,y,gA{x,y)~'^A{y,y)A{x,y)), 

for x,y,y G U and g € G. These equations may be considered as the discrete version 
of the Lagrange-Poincare equations for a G-invariant continuous Lagrangian (see 
[5] for the local expression of the Lagrange-Poincare equations). 

Note that if A : J7 x t/ ^ G is the constant map A{x, y) = e, for all (x, y) G UxU, 
or, in other words, Ad is the trivial connection then equations (|5.9p and ()5.10p may 
be written as 

D2L{{x,y),gk) + DiL{{y, z), gk+i) = 0, ,^ 
lik+i{y,z) = Ad*g^fik{x,y), 

where 

Aife(x,y) = d{r*g^L(^s:^y^ ))(e) 
for {x,y) £U X U (compare equations (|5.1ip with equations (|5.ip and (|5.2p ). 

Discrete Elroy's beanie. As an example of a lagrangian system on an Atiyah grou- 
poid, we consider a discretization of the Elroy's beanie, which is, probably, the most 
simple example of a dynamical system with a non-Abelian Lie group of symmetries. 
The continuous system consists in two planar rigid bodies attached at their centers 
of mass, moving freely in the plane. The configuration space is Q = SE{2) x with 
coordinates {x,y,6,4')i where the three first coordinates describe the position and 
orientation of the center of mass of the first body and the last one the relative ori- 
entation between both bodies. The continuous system is described by a Lagrangian 
Lc{x, y, 0, ^p, X, y, 6, til)) — \m{x^ + y'^) + \li9'^ + ^l2{9 + 'ip)'^ -¥{-4}) where m denotes 
the mass of the system, Ii and I2 are the inertias of the first and the second body, 
respectively, and V is the potential energy. The system admits reduction by SE{2) 
symmetry. In fact, the reduced lagrangian Ic : TQ/SE{2) ~ S"^ x R x se(2) M is 

U^p, ^, ni,n2,n3) = ^minj + nl) + i(/i + i2)nl + Ij^i^' - t^(^) 

Z Z Z li + I2 

where se(2) is the Lie algebra of SE{2), fli = ^1, Q2 — C2, ^^3=^3" TT+h"^ 

and (Ci7C2,'?3) are the coordinates of an element of se(2) with respect to the basis 

/ooi\ /ooo\ , /oio\^^,,,,, . . 

61= 000,62= 001 and 63 = -1 . Note that Si = x cos + y sm 0, 
^Vooo/'^Voooy "^Vooo/ ^ ^ 

h 

^2 = — isin6' + ycosO and ^3 = — — —rp (for more details, see |15 1 127 j ). 

Ii + h 

/cos Ok — sin 6*^ Xk\ 

In order to discretize this system, consider gk — I sinOk cosOk vk € SE{2). We 

V 1 / 

construct the discrete connection Ad : {SE{2) x S^) x {SE{2) x S^) — > SE{2) 
defined by Ad{{gk,ipk), {gk+i,ipk+i)) = .g^+i ^(V'fc, V'fe+Offfc \ where 

/ cos(7j^AV>fc) -sin(7^A,Afc) 0\ 
A{ijk,i'k+l) = sin(7^AV;,) cos(7^A^.) 
\ 1/ 

Here, Aipk = 4'k+i — V'fc- The discrete connection Ad precisely induces the mechan- 
ical connection associated with the iS'£'(2)-invariant metric Q on Q: 

g = mdx ig)dx + mdy ®dy+ [h + I2)d6 ® dO + hdO i® dip + hdip ®d9 + hdip ® dtp 

We remark that the continuous Lagrangian Lc is the kinetic energy associated 
with Q minus the potential energy V . 
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Next, we consider the Atiyah groupoid {Q x Q) / SE{2). As we know, using the 
discrete connection Ad, one may define a local isomorphism between the Atiyah 
groupoid [Q x Q)/SE{2) and the product manifold U x U x SE{2), U being an 
open subset of R. Then, as a local discretization of the reduced Lagrangian Ic, we 
introduce the discrete Lagrangian Id on U x U x SE(2) given by 



1 



ld{lpk,'(pk+l,^{l)k,^{2)k,^{3)k) = T^TI"^ 



2/l2 

1 LL 



^fl)k + ^(2)fe 



^[l-cos(f^(3).)] + 2 7^(^^j -n ^ ) 

where i^ruk ~ Aa;feCos6'fe + AykSiiiOk, ^(2)k — ^Axkshi9k + AykCOsBk and 



fipk + ipk+i 



n 



(3)fc 



-A9k - TifijAVfe 



Now, if we denote by qk = {ipk,iJk+i,^(i)kM{2)k,^{3)k) then 

(0^1 _ = cos(f7(3)fe + ^AV-fe)^ - sin(f7(3), + ^A^,)^ 



(0,CeJ _^ = sin(f](3)fe + T^A^,)^ +cos(l](3)fe + ^AV;,)gj^ 



an, 



an 



9fc 



(2)fc 



/2 



Thus, the reduced Discrete Euler-Lagrange equations 



d 



qk Qk + l Oip qk Oip 



Id^O 



are 



^^(i)fc+i = ^(i)k cos(f](3)fe + j^Aipk) - ^{2)k sin(f^(3)fc + j^Atpk) 
^{2)k+i = f^(i)fe sin(f7(3)fc + j^^Atpk) + ^(2)k cos(f2(3)fe + j^^Aipk) 

^{3)k+l — ^(3)fe 

/i + /2 /i2 ~ 2\di;^ 2 2 

These equations are a discretization of the corresponding reduced equations for the 
continuous system (see |15]). In a forthcoming paper jlO] , we will give a complete 
description of this example comparing with the continuous equations. 



5.5. Reduction of discrete Lagrangian systems. Next, we will present some 
examples of Lie groupoid epimorphisms which allow to do reduction. 

• Let G be a Lie group and consider the pair groupoid G x G over G. Consider 
also G as a groupoid over one point. Then we have that the map 

$i : G X G — > G 

is a Lie groupoid morphism, which is obviously a submersion. Thus, using Corollary 
14. 7|. it follows that the discrete Euler-Lagrange equations for a left invariant discrete 
Lagrangian on G x G reduce to the discrete Lie-Poisson equations on G for the 
reduced Lagrangian. This case appears in [26] as was first noticed by |31| , and also 
appear later in [H H [H HD] . 
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Alternatively, one can do reduction of a right-invariant Lagrangian by using the 
morphism 

^r- GxG — > G 

• Let G be a Lie group acting on a manifold M by the left. We consider a discrete 
Lagrangian on G x G which depends on the variables of M as parameters Lmig, h). 
In general, the Lagrangian will not be invariant under the action of G, that is 
Lm{g, h) ^ Lm{rg, rh). Nevertheless, it can happen that Lm{rg, rh) = L^-i^^g^ h). 
In such cases we can consider the Lie groupoid G x G x M over G x M where 
accordingly one consider the elements in M as parameters. Then the Lagrangian 
can be considered as a function on the groupoid G x G x M given by L{g, h, m) = 
Lm{g,h) so that the above property reads L(rg,rh,rm) — L{g,h,m). Thus we 
define the reduction map 

$ : GxGx M — > GxM 

{g,h,m) {g-^h,g~^m) 

where on G x we consider the transformation Lie groupoid defined by the right 
action m- g — g^^m. Since this map is a submersion, the Euler-Lagrange equations 
onGxGx M reduces to the Euler-Lagrange equations on G x M. This case occurs 
in the Lagrange top that was considered as an example in Section [5T51 (see also [2]). 

• Another interesting case is that of a G-invariant Lagrangian L defined on the pair 
groupoid L : Q x Q — s- M, where p : Q — > M is a G-principal bundle. In this 
case we can reduce to the Atiyah gauge groupoid by means of the map 

QxQ {QxQ)/G 

Thus the discrete Euler-Lagrange equations reduce to the so called discrete Lagrange- 
Poincare equations. 

6. Conclusions and outlook 

In this paper we have elucidated the geometrical framework for discrete Me- 
chanics on Lie groupoids. Using as a main tool the natural Lie algebroid structure 
on the vector bundle tt"^ : V^G —^Gwe have found intrinsic expressions for the 
discrete Euler-Lagrange equations. We introduce the Poincare-Cartan sections, the 
discrete Legendre transformations and the discrete evolution operator in the La- 
grangian and in the Hamiltonian formalism. The notion of regularity has been 
completely characterized and we prove the symplecticity of the discrete evolution 
operators. Moreover, we have studied the symmetries of discrete Lagrangians on 
Lie groupoids relating them with constants of the motion via Noether's Theorems. 
The applicablity of these developments has been stated in several interesting exam- 
ples, in particular for the case of discrete Lagrange-Poincare equations. In fact, the 
general theory of discrete symmetry reduction naturally follows from our results. 

In this paper we have confined ourselves to the geometrical aspects of mechanics 
on Lie groupoids. In a forthcoming paper (see [TU]) we will study the construction 
of geometric integrators for mechanical systems on Lie algebroids. We will intro- 
duce the exact discrete Lagrangian and we will discuss different discretizations of 
a continuous Lagrangian and its numerical implementation. 

Another different aspect we will work on it in the future is to develop natural 
extensions of the above theories for forced systems and systems with holonomic and 
nonholonomic constraints. 
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